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Verfahren zur Erzeugung einer Einhiillenden 



Die Erfindung betrifft ein Verfahren zur automatischen Erzeu- 
gung einer Einhiillenden fur ein rechnerverf ugbares Konstruk- 
tionsmodell, das durch eine vorgegebene Vernetzung in Finite 
Elemente zerlegt ist. 

Unter einer Einhiillenden wird eine rechnerverf iigbare und na- 
herungsweise Beschreibung der Oberflache des vorgegebenen 
Konstruktionsmodells durch eine Menge von Flachenelementen 
verstanden. Die Einhullende approximiert die Oberflache mit 
einer vorgegebenen Maxima labweichung und vernachlassigt L6- 
cher, Aussparungen, Bohrungen, Taschen, Rillen und andere Un~ 
regelm&Sigkeiten, wenn deren Abmessungen kleiner als die vor- 
gegebene Maximalabweichung sind. 

Eine solche Einhiillende wird benotigt, urn physikalische Vor- 
gange an der Oberflache eines technischen Systems durch eine 
Finite-Elemente- Simulation nachzubilden . Hierfur wird das 
technische System durch ein rechnerverf ugbares Konstruktions- 
modell nachgebildet . Beispielsweise ist das modellierte Sys- 
tem ein Motor oder eine Karosserie eines Kraf tf ahrzeugs, und 
die Gerauschemission von der Oberflache des Motors bzw. der 
Karosserie wird durch eine Akustik-Simulation vorhergesagt . 
In der Akustik-Simulation werden nur Vorgange auf der Ober- 
flache und der angrenzenden Umgebung berucksichtigt , nicht 
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aber Vorgange im Inneren des Motors bzw. der Karosserie. Da- 
mit die Rechenzeit der Akustik-Simulation moglichst gering 
ist, approximiert die Einhullende die Oberflache nur so genau 
wie fur die Simulation notig. 

Ein Sof tware-Werkzeug fur Akustik- und andere Simulationen 
namens „ Virtual . Lab" wird unter http://www.lmsintl.com, abge- 
fragt am 24. 10. 2003, beschrieben. Erwahnt wird die Erzeu- 
gung einer Einhullenden, dort „acoustic surface meshing pro- 
cess^ genannt, aus einer vorgegebenen Vernetzung. Aus einer 
Vernetzung eines beispielhaf ten Antriebstrangs mit 150.000 
Finiten Elementen wird eine Einhullende mit 5.000 Flachenele- 
menten erzeugt . Nicht offenbart wird, wie die Einhullende er- 
zeugt wird. 

In DE 10023377 C2 wird ein Verfahren zur Finite-Elemente- 
Simulation beschrieben. Ein vorgegebenes zwei- oder dreidi- 
mensionales Simulationsgebiet wird mit quadratischen Flachen- 
elementen bzw. wiirf elf ormigen Volumenelementen uberdeckt . Er- 
mittelt wird, welche der Gitterzellen ganz, teilweise oder 
gar nicht im Simulationsgebiet liegen. Wie dies geschieht, 
wird nicht offengelegt. Im Simulationsgebiet werden anschlie- 
Send approximierende Funktionen mit Hilfe von B-Splines uber 
den Gitterpunkten definiert, und die Simulation wird mit Hil- 
fe dieser Funktionen durchgefiihrt . 

Aus DE 10056107 Al ist ein Verfahren zur Ermittlung von Klap- 
pergerauschen an Kraf tf ahrzeugen bekannt . Ein Konstruktions- 
modell eines Fahrzeugs wird in Elemente und Gitter unter- 
teilt. Weiterhin wird der Def initionsbereich des Modells in 
Unter-Def initionsbereiche unterteilt, und die Elemente und 
Gitter eines Unter-Def init ionsbereichs werden jeweils separat 
untersucht . 

In DE 19919891 Al wird die Erzeugung eines virtuellen Windka- 
nals z. B. fur ein Kraf tf ahrzeug beschrieben. Ein vorgegebe- 
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nes Konstruktlonsmodell wird vernetzt (eine „Masche" wird er- 
zeugt) , beispxelsweise ein Oberf lachenmodell mit 60.000 drei- 
eckigen Flachenelementen. Fur diese Vernetzung werden vor- 
zugsweise komrnerzielle Werkzeuge verwendet . Bei Bedarf wird 
das Konstruktlonsmodell verfeinert, z. B. in Bereichen mit 
Verwirbelungen . 

Der Erfindung liegt die Aufgabe zugrunde, ein mit geringer 
Rectienzeit durrchf iihrbares Verfahren zur automatischen Erzeu- 
gung einer Einhullenden, welche die Oberflache eines Kon- 
struktionsmodells mit vorgegebener Maximal abweichung approxi- 
miert , berei t zustellen . 

Die Aufgabe wird durch das Verfahren nach Anspruch 1 gelost. 
Vorteilhafte Ausgestaltungen sind in den Unteranspruchen an- 
gegeben. 

Vorgegeben ist eine Vernetzung des Konstruktionsmodells. Die 
Vernetzung besteht aus Finiten Elementen mit Knotenpunkten. 
Das Verfahren umfaJSt folgende Schritte: 

- Ermittelt wird ein Quader, in dem die Vernetzung des Kon- 
struktionsmodells vollstandig enthalten ist. 

- Dieser Quader wird in Volumenelemente zerlegt . Die Zerle- 
gung ist eine vollstandige, so daS jeder Punkt des Quaders 
und damit j eder Punkt des Konstruktionsmodells in ein Vo- 
lumenelement fallt. Der Quader wird so zerlegt, dafi jede 
Kante jedes dieser Volumenelemente kiirzer als oder genauso 
lang wie eine vorgegebene Schranke ist. Diese Schranke ist 
kleiner oder gleich der einzuhaltenden Maximalabweichung 
und kleiner- oder gleich der fur die nachfolgende Simulati- 
on gefordexrte Feinheit der Einhullenden. 

- Fur jedes Volumenelement dieses Quaders wird gepruft, ob 
das Volumenelement mit mindestens einem ermittelten Fini- 
ten Element iiberlappt. Eine Uberlappung liegt bereits vor, 
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wenn das Finite Element und das Volumenelement einen ein- 
zigen Punlct gemeinsam haben. 

- Die Menge der liberlappenden Volumenelemente bildet einen 
geometrischen Korper. Dieser Korper wird ermittelt. 

- Der Korper wird nach auSen durch die auSeren Begrenzungs- 

f lachen der aufienliegenden Volumenelemente begrenzt. Diese 
Begrenzungsf lachen werden ermittelt. 

- Die gesuchte Einhullende wird aus der Menge der dergestalt 
ermittelten Begrenzungsf lachen zusammengesetzt . 

Die erf indungsgemaS erzeugte Einhullende weicht von der zu 
approximierenden Oberflache hochstens urn eine vorgegebene Ma- 
ximal abweichung ab. Die Zerlegung des Quaders, der das Kon- 
struktionsmodell vollstandig enthalt, wird mit einer vorgege- 
benen Genauigkeit durchgef uhrt . Die Zerlegung wird so durch- 
gefiihrt, daS die groSte Kantenlange der Volumenelemente klei- 
ner oder gleich der Maximalabweichung ist. Diese groSte Kan- 
tenlange ist eine obere Grenze fur die Abweichung zwischen 
der Einhullenden, die aus den au&enliegenden Begrenzungsf la- 
chen der uberlappenden Volumenelemente zusammengesetzt wird, 
und der durch die Vernetzung approximierten tatsachlichen O- 
berflache des Konstruktionsmodells . Weil jede Kante kleiner 
oder gleich der vorgegebenen Schranke ist, ist die Maximalab- 
weichung hochstens so groS wie diese vorgegebene Schranke. 
Denn die Kanten der Begrenzungsf lachen sind hochstens so lang 
wie die maximale Kantenlange der Volumenelemente. Weil uber- 
lappende Volumenelemente ausgewahlt werden und keine anderen, 
stellt das Verfahren sicher, daJS der Abstand zwischen der ap- 
proximierenden Vernetzung und der Einhullenden hochstens die 
maximale Kantenlange und damit kleiner oder gleich der vorge- 
gebenen Maximalabweichung ist . 
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Das Verfahren macht es uberfliissig, die vorgegebene Vernet- 
zung des Konstruktionsmodells zu verandern. Dadurch werden 
Rechenschritte und somit Verarbeitungszeit eingespart. Diese 
Einsparung ist urn so grower, je feiner die vorgegebene Ver- 
netzung im Vergleich zur Maximal abweichung ist. Das Verfahren 
erfordert bei N Quadern und M Finiten Elementen maximal M * N 
Priifschritte. 

Das erf indungsgemafie Verfahren hangt nicht von der Art ab, 
wie das Konstruktionsmodell vernetzt wurde. Das Verfahren 
laSt sich insbesondere sowohl dann anwenden, wenn die Vernet- 
zung des Konstruktionsmodells aus Flachenelementen besteht, 
als auch dann, wenn sie aus Volumenelementen besteht . Das 
Verfahren laSt sich auch fur jede Form der Flachenelemente 
anwenden. Der Quader wird vollig unabhangig von der Vernet- 
zung des Konstruktionsmodells in Volumenelementen zerlegt. 
Dadurch wird die Einhullende so erzeugt, wie die nachfolgende 
Simulation dies erfordert. Diese nachfolgende Simulation kann 
eine andere sein als diejenige; fur welche die vorgegebene 
Vernetzung erzeugt wurde . Beispielsweise wurde die Vernetzung 
fur eine Simulation der im System auftretenden Spannungen er- 
zeugt, die Einhullende wird hingegen fur eine Akustik- 
Simulation verwendet. 

Nur diejenigen Begrenzungsf lachen der uberlappenden Volumen- 
elemente werden ausgewahlt, die den durch die uberlappenden 
Volumenelemente gebildeten Korper nach aufien begrenzen. Da- 
durch wird erreicht, daS die innere Struktur des Konstrukti- 
onsmodells bei der Erzeugung der Einhullenden nicht beruck- 
sichtigt wird. Weiterhin wird erreicht, daS die erf indungsge- 
ma£ erzeugte Einhullende ist daher immer geschlossen und re- 
gelmafiig ist. 

Falls das vorgegebene Konstruktionsmodell Locher aufweist, z. 
B. Bohrungen oder Taschen, so hangt es von deren Abmessungen 
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ab, ob Approximationen dieser Locher in der Einhullenden auf- 
treten oder ob die Locher bei der Approximation nicht beruck- 
sichtigt werden. Falls der Durchmesser des Lochs groSer ist 
als die gro&te Kantenlange der Volumenelemente des Quaders, 
so liegt mindestens ein Volumenelement des Quaders ganz oder 
teilweise in dem Loch, und Begrenzungsf lichen dieses Quaders 
sind auSenliegende und gehoren zur Einhullenden. Ansonsten 
wird das Loch vernachlassigt. Durch die frei wahlbare Grofie 
der Volumenelemente des Quaders wird also gesteuert, welche 
Locher beriicksichtigt und welche vernachlassigt werden. 

Die erzeugte Einhullende besitzt wesentlich weniger Knoten- 
punkte als die vorgegebene Vernetzung des Konstruktionsmo- 
dells. Eine Finite-Elemente-Simulation mit den Knotenpunkten 
der Einhullenden verwendet ein Gleichungssystem mit weniger 
Unbekannten und erfordert daher weniger Rechenzeit. Diese 
Einsparung fallt insbesondere dann ins Gewicht, wenn die Si- 
mulation mehrmals durchgefuhrt wird, z. B. urn verschiedene 
Konstruktionsstande miteinander zu vergleichen. 

Die Ausgestaltung nach Anspruch 2 erfordert besonders wenig 
Rechenzeit. Der Quader wird so zerlegt, daS die kiirzeste Kan- 
te jedes Volumenelement s groSer oder gleich der langsten Kan- 
te jedes Finiten Elements der vorgegebenen Vernetzung ist. 
Die Zerlegung des Quaders ist frei wahlbar, und die Genauig- 
keit, mit der die Einhullende die Oberflache des Konstrukti- 
onsmodells approximieren soli, wird durch die Anf orderungen 
der Simulation, fur welche die Einhullende verwendet wird, 
vorgegeben. Daher ist es meistens moglich, daS die Volumen- 
elemente des Quaders langere Kanten als die Flachenelemente 
der Vernetzung des Konstruktionsmodells haben. 

Die Priifung, ob ein Finites Element der Vernetzung und ein 
Volumenelement des Quaders uberlappen oder nicht, wird durch 
die Ausgestaltung nach Anspruch 2 auf besonders einfache Wei- 
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se durchgef uhrt . Gepruft wird, ob mindestens ein Knotenpunkt 
des Finiten Elements im Volumenelement liegt. Wird ein sol- 
cher Knotenpunkt gefunden, so wird die Prufung abgebrochen, 
und entschieden wird, daS das Volumenelement mit dem Finiten 
Element uberlappt . Falls kein Knotenpunkt des Finiten Ele- 
ments im Volumenelement liegt, so kann auch kein anderer 
Punkt im Volumenelement liegen. Dies folgt aus der Ausgestal- 
tung, dafi die Kanten des Volumenelement s langer als die des 
Finiten Elements sind. Weil kein anderer Punkt des Finiten 
Elements im Volumenelement liegt, liberlappen Finites Element 
und Volumenelement nicht. Weil nur die Knotenpunkte unter- 
sucht werden und keine anderen Punkt e, wird besonders wenig 
Rechenzeit benotigt . 

Im Folgenden wird ein. Ausf uhrungsbeispiel der Erfindung naher 
beschrieben, dabei zeigen: 

Fig. 1. Knotenpunkte einer vorgegebenen unregelmaSigen Fla- 
che; 

Fig. 2. ein die Flache enthaltendes Rechteck; 

Fig. 3. eine Zerlegung des Rechtecks in kleine Rechtecke; 

Fig. 4. die Rechtecke der Zerlegung, die Knotenpunkte enthal- 
ten; 

Fig. 5. die entstehende massive zweidimensionale Flache; 
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Fig. 6. die erzeugte Einhullende. 

Das erf indungsgemaSe Verfahren wird zunachst vereinfacht an 
einem zweidimensionalen Beispiel veranschaulicht . In diesem 
Beispiel ist eine Einhullende fur eine unregelmaSige Flache 
in der Ebene zu erzeugen. Diese unregelmaSige Flache tritt in 
dieser Veranschaulichung an die Stelle des vorgegebenen Kon- 
struktionsmodells . Die zu erzeugende Einhiillende ist ein 
Streckenzug in dieser Ebene. Fig. 1 zeigt einige der Knoten- 
punkte derjenigen Flachenelemente , die zu einer Vernetzung 
der Flache gehoren. Der Ubersicht lichkeit halber sind viele 
Knotenpunkte nicht dargestellt. Die durch die Knotenpunkte 
gebildeten Flachenelemente der unregelmaSigen Flache sind e- 
benfalls nicht dargestellt. 

Ein Rechteck wird ermittelt, in dem die vorgegebene unregel- 
maSige Flache und damit die Vernetzung enthalten ist. Die un- 
regelmaSige Flache und das ermittelte Rechteck werden in 
Fig. 2 gezeigt. Dieses Rechteck tritt an die Stelle des Qua- 
ders . 

Das Rechteck wird in glei'ch groSe Rechtecke mit vorgegebenen 
Kantenlangen zerlegt. Fig. 3 zeigt diese Rechtecke der Zerle- 
gung, die an die Stelle der Volumenelemente treten, sowie die 
Knotenpunkte der Fig. 1. 

Gepriift wird, welche Rechtecke der Zerlegung jeweils mindes- 
tens einen Knotenpunkt der unregelmaSigen Flache umfassen. 
Diese uberlappenden Rechtecke sind in Fig. 4 schraffiert dar- 
gestellt . 

Die in Fig. 4 schraffiert dargest ellten Rechtecke bilden eine 
zweidimensionale Flache. Die Hohlraume in dieser zweidimensi- 
onalen Flache werden ausgefullt . Dadurch entsteht die in 
Fig. 5 gezeigte massive zweidimensionale Flache. 
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Die Begrenzungslinien dieser massiven Flache werden zur ge- 
suchten Einhullenden zusammengesetzt . Diese Begrenzungsli- 
nien, die einen Streckenzug bilden, werden in Fig. 6 durch 
fette Linien gezeigt. 

Das nun folgende Ausf iihrungsbei spiel bezieht sich auf einen 
Bestandteil eines neuen Kraf tf ahrzeuges als dem technischen 
System. Die erzeugte Einhiillende wird beispielsweise dazu 
verwendet, die Schallemission von einem festen Kdrper, z. B. 
einem Motor oder einer Karosserie eines Kraf tf ahrzeuges , 
durch eine Finite-Elemente- Simulation zu untersuchen. Die 
Schallemission wird durch Korperschall hervorgeruf en . Dieser 
entsteht durch Zug- und/oder Schubspannungen in Festkorpern. 

Die Schallemission wird durch Vibrationen an der Oberflache 
des Motors oder der Karosserie hervorgeruf en . Diese Vibratio- 
nen wirken auf die umgebende Luft ein. Die Vibrationen werden 
in einen Anteil senkrecht zur Oberflache und einen parallel 
zur Oberflache zerlegt. Nur die Anteile der Vibrationen senk- 
recht zur Oberflache werden in der Akustik-Simulation beruck- 
sichtigt, nicht aber die Anteile parallel zur Oberflache. 

Aus der Strukturdynamik des Motors und seiner Oberflache re- 
sultiert eine Verteilung der Schnelle der Vibrationen. Diese 
Schnelle wird durch eine Akustik-Simulation ermittelt . Mit 
der Schnelle wird die Auslenkungsgeschwindigkeit ds/dt be- 
zeichnet, wobei s(t) die Auslenkung der Schallwelle senkrecht 
zur Oberflache in Abhangigkeit von der Zeit t ist und wobei 
gilt: 

ds/dt = -co * A * sin (cot + kx + . 

Die Strukturdynamik des Motors wird mit einer ersten Finite - 
Elemente-Simulation berechnet . Die berechnete Strukturdynamik 
umfaSt die Verschiebung, die Geschwindigkeit oder die Be- 
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schleunigung z. B. der Verschiebung der Oberflache. Aus die- 
ser wird die Schallschnelle der Oberflache berechnet . 

Vorgegeben ist ein rechnerverf ugbares Konstruktionsmodell ei- 
nes zu konstruierenden oder zu untersuchenden Systems, z. B. 
des Motors. Dieses Konstruktionsmodell ist auf einem Daten- 
speicher einer Datenverarbeitungsanlage abgespeichert und hat 
vorzugsweise die Form eirtes dreidimensionaleri CAD-Modells. 
Das Konstruktionsmodell wird beispielsweise mit Hilfe von 
Punkten, Kurven, Vektoren, Linien, Konturen, polygoniale Net- 
ze, gekrummten Oberflachen und/oder Volumenelementen in dem 
Datenspeicher reprasentiert . Das erf indungsgemaSe Verfahren 
wird mit dieser. Datenverarbeitungsanlage oder einer anderen, 
die Lesezugriff auf den Datenspeicher hat, durchgef uhrt . 

Die Methode der Finiten Elemente ist aus „Dubbel - Taschen- 
buch fur den Maschinenbau" , 20. Auflage, Springer- Verlag, 
2001, C 48 bis C 50, aus B. Klein: „FEM - Grundlagen und An- 
wendungen der Finite-Elemente-Methode" , Vieweg-Verlag, 
3. Auflage, 1999, aus T. R. Chandrupalta & A. D. Belegundu: 
„ Introduction to Finite Element in Engineering", Prentice- 
Hall, 1991, sowie aus DE 19927941 CI bekannt . Durch eine Fi- 
nite-Elemente-Simulation werden physikal ische Vorgange unter- 
schiedlicher Arten untersucht, beispielsweise Schwingungs vor- 
gange, z. B. die Gerauschemission von der Oberflache eines 
Motors oder einer Karosserie eines Kraft fahrzeuges . 

In dem vorgegebenen Konstruktionsmodell wird eine bestimmte 
Menge von Punkten festgelegt, die Knotenpunkte heiSen. Als 
Finite Elemente werden diejenigen Flachen- oder Volumenele- 
mente bezeichnet, deren Ecken durch Knotenpunkte definiert 
werden. Die Knotenpunkte bilden ein Netz in dem Konstrukti- 
onsmodell, weswegen der Vorgang, Knotenpunkte festzulegen und 
Finite Elemente zu erzeugen, Vernetzen und die Finiten Ele- 
mente Vernetzung des Modells genannt werden. In DE 
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10010408 Al und DE 19933314 Al werden Verfahren zum Verne tzen 
eines gegebenen Konstruktionsmodells beschrieben. 

Moglich ist, das Konstruktionsmo'dell beim Vernetzen vollstan- 
dig in Finite Elemente zu zerlegen. In diesem Fall gehort je- 
der Punkt des Konstruktionsmodells zu mindestens einem Fini- 
ten Element. Oft reicht aber eine Naheruragslosung aus ; die 
Knotenpunkte einspart und damit zu einem kleineren Glei- 
chungssystem f uhrt . Falls beispielsweise die Karosserie eines 
Kraftf ahrzeugs untersucht wird, so reicht es meist aus, die 
Blechteile der Karosserie durch Flachenelemente naherungswei- 
se zu vernetzen. Fiigeverbindungen zwischen Blechteilen, z. B. 
SchweiS- oder Klebenahte, werden mitt els Volumenelementen 
vernetzt. Dadurch lassen sich Vorgange in. den Fiigeverbindun- 
gen, z. B. Verschiebungen und Querspannungen, genauer unter- 
suchen. In manchen Anwendungen reicht es aus, die Oberflache 
des Konstruktionsmodells zu vernetzen, beispielsweise durch 
eine „Triangulierung w mittels dreieckigerr Flachenelemente. 
Das erf indungsgemaSe Verfahren laSt sich ohne Modifikation 
fur jede dieser Vernetzungen anwenden. Eine beliebige Vernet- 
zung des Konstruktionsmodells wird dem Verfahren vorgegeben. 

Je nach Problemstellung werden die Verschiebungen dieser Kno- 
tenpunkte und/oder Rotationen der Finiten Elemente in diesen 
Knotenpunkten und/oder die Spannungen in bestimmten Punkten 
dieser Finiten Elemente, namlich in den Integrationspunkten, 
als Unbekannte eingefuhrt. Gleichungen werden aufgestellt, 
welche z. B. die Verschiebungen, Rotationen und/oder Spannun- 
gen innerhalb eines Finiten Elements nahe rungs we ise beschrei- 
ben. Weitere Gleichungen result ieren aus Abhangigkeiten zwi- 
schen verschiedenen Finiten Elementen. Insgesamt wird ein oft 
sehr umf angreiches Gleichungssystem mit d.en Knotenpunkt- 
Verschiebungen, Knot enpunkt -Rotationen, Spannungen in Integ- 
rationspunkten oder weitere GroSen als Unbekannte aufgestellt 
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und numerisch geldst. Aus der Losung des Gleichungssystems 
lassen sich Aussagen uber das physikalische Verhalten des 
Systems gewinnen. 

Die Einhullende laSt sich auch fur eine Simulation gemafi der 
Boundary-Elemente-Methode verwenden. Diese Methode wird in 
Dubbel, a.a.O., C5 0 - C52, beschrieben. Wie eine Akustik- 
Simulation mit der Finite-Elemente-Methode oder der Boundary- 
Elemente-Methode durchgefuhrt wird, ist z. B. in Dubbel, 
a.a.O. , 033 - 035, beschrieben. 

Zur Durchfuhrung des erf indungsgemaSen Verfahrens wird ledig- 
lich die Vernetzung des Konstruktionsmodells benotigt. Weite- 
re Inf ormationen uber das Konstruktionsmodell , z. B. Kriim- 
mungsradien oder Fasen oder die Oberf lachenbeschaf f enheit , 
erfordert das Verf ahren .nicht . Die Datenverarbeitung, die das 
Verfahren durchfiihrt, liest die Vernetzung des Konstruktions- 
modells ein. Hierbei werden die im Folgenden beschriebenen 
Inf ormationen liber Knotenpunkte und uber die Finiten Element e 
eingelesen. 

Ein Knotenpunkt ist durch eine eindeutige Kennung, z. B. eine 
fortlaufende Nummer, und seiner Position, gekennzeichnet . Die 
Kennung ist in der Vernetzung eindeutig, weswegen zwei ver- 
schiedene Knotenpunkte stets zwei verschiedene Kermungen ha- 
ben. Die Position des Knotenpunkts wird durch je eine x- , y- 
und z-Koordinate in einem dreidimensionalen Koordinatensystem 
gekennzeichnet. Beispielsweise urn Verschiebungen und Verfor- 
mungen unter dem EinfluS eines Lastfalls zu untersuchen, kann 
vorgesehen sein, daJS zwei Knotenpunkte zweier verschiedener 
Finiten Elemente vor der Verformung dieselbe Position und - 
daher dieselben x-, y- und z-Koordinaten haben. Diese beiden 
Knotenpunkte werden durch verschiedene Kermungen unterschie- 
den. 
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Man unterscheidet zweidimensionale Finite Elemente, die im 
folgenden Flachenelemente genannt werden, und dreidimensiona- 
le Finite Elemente, die im Folgenden Volumenelemente hieiSen. 
Gekrummte Flachen und Korper, die durch gekrummte Flachen im 
Raum angenahert werden, werden oft mittels Flacheneletnenten 
vernetzt. Die gebrauchlichsten Flachenelemente sind Dr-eiecke 
und Vierecke. Alle Ecken der Flachenelemente sind Knotenpunk- 
te. Manchmal werden Vierecke mit sechs Knotenpunkten verwen- 
det. Neben rechteckigen Flachenelementen konnen auch andere 
viereckige Flachenelemente vorgesehen sein. Flachenelemente 
werden vorzugsweise durch eine isoparametrische Darstellung 
beschrieben, die z. B. aus T. R. Chandrupalta & A. D. Bele- 
gundu, a.a.O., bekannt ist. 

Als Volumenelemente werden in der Regel Wurfel und anciere 
Quader verwendet, also Volumenelemente mit sechs ausschliefc- 
lich rechteckigen Begrenzungsf lachen. Derartige Volumenele- 
mente werden oft als Hexaeder bezeichnet. TypischerwedLse sind 
die acht Ecken des Quaders Knotenpunkte des Volumenelements. 
Verwendet werden Volumenelemente mit acht bis zwanzig Knoten- 
punkten. 

Jedes Finite Element ist durch eine in der Vemetzung eindeu- 
tige Kennung, z. B. eine fortlaufende Nummer, und durch seine 
Knotenpunkte gekennzeichnet . Vorzugsweise umfassen die Infor- 
mationen uber das Finite Element eine Auflistung der Kennun- 
gen seiner Knotenpunkte sowie eine Kennzeichnung, von welcher 
Art das Finite Element ist. 

Als ersten Schritt bei der Durchfuhrung des Ausfuhrungsbei- 
spiels liest die Datenverarbeitung, die das Verfahren aus- 
fiihrt, die Vernetzung ein. Ein Quader wird ermittelt , in dem 
die Vernetzung des Konstruktionsmodells vollstandig erathalten 
ist. Vorzugsweise wird diese Quader so gelegt, dafi jec3.e sei- 
ner Begrenzungsf lachen senkrecht auf einer Koordinatenachse 
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liegt und somit seine Kanten achsenparallel ausgerichtet 
sind. Jede der sechs Begrenzungsf lachen liegt dann in einer 
der folgenden Ebenen: 

{ (x,y,z) | x m x_min }, { (x,y,z) | x = xjnax }, 
{ (x,y,z) | y = y_min }, { (x,y,z) | y = yjnax } # 
{ (x,y,z) | z = z_min }, { (x f y,z) | z - z_max }. 

Der Wert x_min wird mit Hilfe der Positionen der Knotenpunkte 
der eingelesenen Vernetzung berechnet, namlich. als der 
kleinste Wert aller x-Koordinaten der Knotenpunkte. Entspre- 
chend wird x_max als der grofite Wert aller x-Koordinaten der 
Knotenpunkte berechnet. Auf analoge Weise werden y_min, 
y_jnax, z_min und z_max berechnet. Alle Finiten Elemente der 
vorgegebenen Vernetzung liegen innerhalb dieses Quaders, weil 
jeder Eckpunkt eines Finiten Elements ein Knotenpunkt der 
Vernetzung ist und weil die Berechnung von x_min, ... , z_max 
sicherstellt, daS alle Knotenpunkte innerhalb des Quaders 
liegen. 

Der Quader, in dem die Vernetzung enthalten ist, wird in Vo- 
lumenelemente zerlegt. Vorzugsweise haben alle Volumenelemen- 
te die gleichen Abmessungen und die gleiche geometrische Form 
und grenzen nahtlos aneinander. Beispielsweise wird der Qua- 
der in gleich groSe Wurfel zerlegt. In dem Falle, daS ein 
harmonischer Schwingungsvorgang durch eine Finite- Elemente - 
Simulation auf Basis der Einhullenden untersucht werden soli, 
wird wie folgt festgelegt, wie der Quader in Wurfel zerlegt 
wird: Eine halbe Welle des Schwingungsvorgangs wird in der 
Finite-Elemente-Simulation durch einen Streckenzug mit drei 
gleich langen Strecken approximiert , eine voile Welle demnach 
durch sechs Strecken. Daher wird der Quader in Wurfel zer- 
legt, deren Kantenlange ein Sechstel der Wellenlange, also 
der Lange der vollen Welle, ist. Diese Wellenlange hangt von 
der Frequenz des Schwingungsvorgangs und von dem Medium, in 
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dem die Schwingung sich ausbreitet, ab, und ist daher be- 
kannt . 

Gepruft wird, ob die kiirzeste Kantenlange der Volumenelemente 
des Quaders groSer ist als die langste Kante der Vernetzung 
des Konstruktionsmodells. Ist dies nicht der Fall, so wird 
die Vernetzung des Konstruktionsmodells verfeinert. 

Vorzugsweise werden zunachst die Kantenlangen der Volumenele- 
mente in x-, y- und z-Richtung ermittelt oder vorgegeben. An- 
schlieiSend werden die Langen der drei Kanten des Quaders in 
x-, y- und z-Richtung durch die vorgegebenen oder ermittelfcen 
Kantenlangen der Volumenelemente dividiert. Die drei Quotien- 
ten werden auf die nachsten naturlichen Zahlen N_x, N_y unci 
N_z aufgerundet. Dadurch wird die Anzahl der Volumenelemente 
in x-, y- und z-Richtung berechnet . Eine tatsachliche Kanten- 
lange eines Volumenelements kann daher von der vorgegebenen 
leicht abweichen. 

Die Volumenelemente, in die der Quader zerlegt wird, werden 
fortlaufend durchnumeriert . Dadurch erhalt jedes dieser Volu- 
menelemente eine eindeutige Kennung. Vorzugsweise ist die 
Kennung ein Tripel (n_x, n_y, n_z) , das die relative Position 
des Volumenelements im Quader angibt . Das Volumenelement , das 
den Eckpunkt (x_min, yjmin, z_ min) des Quaders enthalt, er- 
halt die Kennung (1, 1, 1) . Das Volumenelement , das den Eck- 
punkt (xjnax, y_max, z_max) des Quaders enthalt, erhalt die 
Kennung (N_x, N_y, N_z) , wobei N_x, N_y und N_z die Anzahl 
der Volumenelemente in x- , y- bzw. z-Richtung angibt. 

Als nachstes wird ermittelt, welche Volumenelemente des Qua- 
ders mit Finiten Elementen der Vernetzung uberlappen. Hierzu 
wird eine Datenstruktur in Form einer dreispaltigen Tabelle 
erzeugt . Anstelle einer Tabelle kann auch eine verkettete 
Liste als Datenstruktur verwendet werden, wobei jedes Listen- 
element drei Eintrage ermoglicht. 
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In der ersten Spalte der Tabelle werden die Kennungen der Vo- 
lumenelemente des Quaders eingetragen, in die zweite Spalte 
in jeder Zeile eine Null. AnschlieSend wird eine Datenstruk- 
tur, in der die Knotenpunkte der Vernetzung abgespeichert 
sind, durchlaufen. Fur jeden Knotenpunkt der Vernetzung wird 
ermittelt, in welches Volumenelement des Quaders diesear Kno- 
tenpunkt fullt. Wie oben dargelegt, liegt jeder Knotenpunkt 
in dem Quader, und der Quader wird vollstandig und ohne Rest 
in Volumenelemente zerlegt . Moglich ist, daS ein Knotenpunkt 
auf der Begrenzungsf lache zweier aneinander angrenzenden Vo- 
lumenelement en liegt. In der zweispaltigen Tabelle wird fur 
die Volumenelemente, in die der Knotenpunkt fallt, der alte 
Eintrag in der zweiten Spalte um Eins vergroSert . Nach Ende - 
des Schrittes gibt die Tabelle fur jedes Volumenelement an, 
wie viele Knotenpunkte in ihm liegen. 

Die Menge der iiberlappenden Volumenelemente ist genau die 
Menge derjenigen Volumenelemente, die mindestens einen Kno- 
tenpunkt enthalten und in der zweiten Spalte der Tabelle ei- 
nen Wert groSer Null besitzen. Diese Menge bildet einen. drei- 
dimensionalen Korper, der vollstandig im Quader liegt. Dieser 
Korper kann Hohlraume enthalten, die im nachsten Schritt ge- 
fiillt werden. Nach Auffullung der Hohlraume entsteht ein mas- 
siver Korper. Die dritte Spalte der Tabelle legt fest, welche 
Volumenelemente des Quaders zu diesem massiven Korper gehoren 
und welche nicht. Hierzu wird in die dritte Spalte der Tabel- 
le eine Eins eingetragen, wenn in der zweiten Spalte ein Wert 
grofier Null steht . 

Die Volumenelemente, die keinen Knotenpunkt enthalten u.nd da- 
her in der zweiten Spalte eine Null haben, gehoren entweder 
zu Hohlraumen im Korper gehoren oder liegen auSerhalb des 
Korpers. Um diese beiden Falle zu unterscheiden, wird wie 
folgt vorgegangen: Fur jedes Volumenelement VE, das keinen 
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Knotenpunkt enthalt und fur das in der dritten Spalte noch 
kein Eintrag vorgenommen wurde, werden die sechs angrenzenden 
Volumenelemente ermittelt, also diejenigen, die mit einer Be- 
grenzungsf lache an VE angrenzen. Ermittelt wird, welche die- 
ser sechs angrenzenden Volumenelemente mindestens einen Kno- 
tenpunkt enthalten und welche keinen. Fur die an VE angren- 
zenden Volumenelemente, die keinen Knotenpunkt enthalten, 
werden wiederum die angrenzenden Volumenelemente ohne Knoten- 
punkte ermittelt, die nicht an VE angrenzen. Diese Ermittlun- 
gen werden fortgesetzt, bis entweder ein zusammenhangender 
Hohlraum, der VE enthalt, ermittelt wurde oder bis der Rand 
des Quaders erreicht ist. Alle Volumenelemente eines zusam- 
menhangenden Hohlraums werden dadurch „gefullt", daS in der 
dritten Spalte jeweils eine Eins eingetragen wird. Wird hin- 
gegen ausgehend von VE der Rand des Quaders erreicht, so 
liegt VE nicht im Inneren des Korpers, und in der dritten 
Spalte werden fur VE und fur die angrenzenden Volumenelemente 
ohne Knotenpunkte je eine Null eingetragen. 

Das Auffullen von Hohlraumen laSt sich auch auf anderem Wege 
durchfuhren. Zunachst wird in der dritten Spalte fur alle Vo- 
lumenelemente mit je mindestens einem Knotenpunkt je eine 
Eins eingetragen. AnschlieSend werden nacheinander alle Volu- 
menelemente ermittelt, die am Rande des Quaders liegen, d. h. 
die mindestens eine Begrenzungsf lache besitzen, die in einer 
Begrenzungsf lache des Quaders liegt. Diejenigen dieser Volu- 
menelemente, die keinen Knotenpunkt der Vernetzung enthalten, 
sind aufienliegende Volumenelemente und erhalten eine Null in 
der dritten Spalte. Diejenigen Volumenelemente, die an diese 
Volumenelemente am Rande angrenzen und ebenfalls keine Kno- 
tenpunkte enthalten, werden ermittelt und erhalten ebenfalls 
eine Null in der dritten Spalte. Wiederholt werden angrenzen- 
de Volumenelemente ohne Knotenpunkte ermittelt, bis alle Vo- 
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lumenelemente untersucht sind und die dritte Spalte vollstan- 
dig ausgef ullt ist . 

Nach dem Auffullen liegt eine aus Volumenelementen des Qua- 
ders bestehende Vernetzung eines massiven dreidimens±onalen 
Korpers vor. Diejenigen Begrenzungsf lachen dieser Vernetzung 
werden ermittelt, die diesen Korper nach auSen begrenzen. 
Vorzugsweise wird hierzu aus dem Volumenmodell , das ist die 
Vernetzung des massiven Korpers mit Volumenelementen, ein O- 
berf lachenmodell, das ist eine Vernetzung der Oberf lsiche mit 
Flachenelementen, erzeugt . Dieser Schritt wird z. B. mit ei- 
nem Praprozessor fur Finite-Elemente-Werkzeuge durchgefxihrt , 
Als erstes wird eine Beschreibung der aus Volumenelementen 
bestehenden Vernetzung des Korpers in einer Datei abgespei- 
chert, z. B. im Datenformat „NASTRAN Bulk Data". Der Prapro- 
zessor „ MEDINA" vermag eine Vernetzung aus einer Datei einzu- 
lesen, wenn diese in einem standardisierten Datenformat, z. 
B. in NASTRAN Bulk Data, vorliegt. Eine Beschreibung von 
MEDINA ist unter 

http : //www. c3pdm. com/des/products/medina/documentation/medina 
-DINA4_e.pdf, abgefragt am 27. 11. 2003, verfugbar. MEDINA 
besitzt eine Funktionalitat, urn aus einem Volumenmodell ein 
Oberf lachenmodell zu erzeugen. Dieses Oberf lachenmodell ist 
die gesuchte Einhullende. 

Vorzugsweise wird das Oberf lachenmodell geglattet („orela- 
xiert") , wobei Treppenstuf en in der Vernetzung entfearnt wer- 
den. Auch die Glattung wird z. B. mit einer entsprechenden 
Funktionalitat von MEDINA durchgef uhrt . 

Die erf indungsgemafi erzeugte Einhullende lalSt sich aiach dazu 
verwenden, Ergebnisse der Finite-Elemente-Simulation, z. B. 
das Schwingungsverhalten oder die Schallausbreitung , zu visu- 
alisieren. Fur diese Visual is ierung werden nur die berechne- 
ten physikalischen GroSen an der Oberflache des untearsuchten 
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Systems benotigt, nicht aber die innere Struktur. Mit cier 
Einhullenden wird eine Datenreduktion durchgef uhrt . Die Be- 
rechnungsergebnisse werden auf die Flachenelemente der Ein- 
hullenden umgerechnet , z. B. indem per Interpolation die Wer- 
te der GroSen in den Eckpunkten dieser Flachenelemente be- 
rechnet werden. Die berechneten GroSen werden vorzugsweise 
auf einem Bildschirm einer Datenverarbeitungsanlage oder auf 
einem Papierausdruck visualisiert . Infolge der Datenreduktion 
laSt sich die Visual isierung schneller erzeugen und veran- 
dern, z. B. wenn ein Benutzer die Visual isierung auf dem 
Bildschirm drehen oder vergroSern will- 
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Patent anspriiche 



1. Verfahren zur automatischen Erzeugung einer Einhullenden, 
welche die Oberflache eines rechner*verf ugbaren Konstruk- 
tionsmodells eines technischen Systems approximiert , 

wobei eine aus Finiten Element en mit Knotenpunkten beste- 
hende Vernetzung des Konstruktionsmodells vorgegeben ist 

und das Verfahren die Schritte umf a£t : 

- Ermitteln eines Quaders, in dem die Vernetzung voll- 
standig enthalten ist, 

- Zerlegen des Quaders in Volumenelemente mit Kanten, 
die hochstens genauso lang^wie eine vorgegebene 
Schranke sind, 

- Fur jedes Volumenelement Prufen, ob das Volumenelement 
mit mindestens einem Finiten Element der Vernetzung 
uberlappt , 

- Ermitteln eines Korpers, der durch die uberlappenden 
Volumenelemente gebildet wird, 

- Ermitteln derjenigen Begrenzungsf lachen, die den Kor- 
per nach aufien begrenzen, und 

- Zusammensetzen der Einhullenden aus den ermittelten 
Begrenzungsf lachen . 
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2. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet , 
daS 

- der Quader so zerlegt wird, daS die kurzeste Kante je- 
des Volumenelements grower oder gleich der langsten 
Kante jedes Finiten Elements ist 

- dann, wenn es einen Knotenpunkt eines Finiten Elements 
gibt, der in einem Volumenelement liegt, Entscheiden, 
da£ das Volumenelement mit dem Finiten Element iiber- 
lappt , und 

- dann, wenn kein Knotenpunkt eines Finiten Elements in 
einem Volumenelement liegt, Entscheiden, daS dieses 
Volumenelement mit dem Finiten Element nicht uber- 
lappt . 

3 . Verfahren nach Anspruch 1 oder Anspruch 2 , 
dadurch gekennzeichnet, dalS 

- die Finiten Elemente der Vernetzung Rechtecke 

- und die Volumenelemente * Quader sind. 

4. Verfahren nach einem der Anspruche 1 bis 3, 
dadurch gekennzeichnet, 
daS 

- eine Wellenlange eines harmonischen Schwingungsvor- 
gangs vorgegeben ist, 
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- als Volumenelemente Wurfel verwendet werden, deren 
Kantenlange ein vorgegebener Bruchteil der Wellenlange 
ist, und 

- unter Verwendung der Einhullenden eine Finite- 
Elemente-Simulation des Schwingungsvorgangs durchge- 
fuhrt wird. 

5. Verfahren nach einem der Anspriiche 1 bis 4, 
dadurch gekennzeichnet , 

daS beim Ermitteln der Begrenzungsf lachen 

- diejenigen Volumenelemente des Quaders ermittelt wer- 
den, die Hohlraume im Inneren des ermittelten Korpers 
bilden, 

- diese Hohlraume gefullt werden und 

- die Begrenzungsf lachen des dadurch erzeugten massiven 
dreidimensionalen Korpers ermittelt werden. 

6. Computerprogramm-Produkt , das direkt in den internen 
Speicher eines Computers geladen werden kann und Soft- 
wareabschnitte umfafit, mit denen ein Verfahren nach einem 
der Anspriiche 1 bis 5 ausgefuhrt werden kann, wenn das 
Produkt auf einem Computer lauft. 

7. Computerprogramm-Produkt, das auf einem von einem Compu- 
ter lesbaren Medium gespeichert ist und das von einem 
Computer lesbare Programm-Mittel aufweist, die den Compu- 
ter veranlassen, ein Verfahren nach einem der Anspriiche 1 
bis 5 auszufuhren. 
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